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(...) 

Rolling lubrication with spherical molecules working as "nanobearings" has failed experimen- 
tally so far, without a full understanding of the physics involved and of the reasons why. Past 
model simulations and common sense have shown that molecules can only roll when they are 
not too closely packed to jam. The same type of model simulations now shows in addition that 
molecular rolling friction can develop deep minima once the molecule's peripheral "pitch" can 
match the substrate periodicity, much as ordinary cogwheels do in a rack-and-pinion system. 
When the pinion-rack matching is bad, the driven molecular rolling becomes discontinuous 
and noisy, whence energy is dissipated and friction is large. This suggests experiments to 
be conducted by varying the rack-and-pinion matching. That could be pursued not only by 
changing molecules and substrates, but also by applying different sliding directions within the 
same system, or by applying pressure, to change the effective matching. 

Keywords: nanotribology; rolling friction; fullerenes 



1. Introduction 

The problem of designing nanomechanical devices, in particular, of reducing fric- 
tion by means of nano- and micro-bearings [l|, HJ is nowadays a hot one. The idea 
of using spherical molecules such as fullerene Cgo as molecular nanobearings gave 
rise to experimental attempts 0, 0, El as well as Molecular Dynamics (MD) simu- 
lations [1, 0, @, [§] • Other potentially rolling molecules were considered as additives 
in oil lubricants, and predicted to provide interesting tribological properties [Io| . 
The simulations indicated that ball-shaped molecules may either slide or rotate 
over a surface, depending on the substrate and the position of the molecule, and 
suggested that the rolling configuration should indeed attain extremely low fric- 
tion 0, \§L with a predicted friction coefficient of the order /i ~ 0.01 — 0.02 [?J or 
smaller [8(. All experimental attempts however met with scarce success so far. A 
single Ceo molecule confined between solid substrates should begin to roll under a 
torque of order 10~ 19 Nm (ll| . However, Cgo molecules generally condense to form 
close-packed islands, their two dimensional (2D) crystalline structure exhibiting 
order at low temperatures. Through a first-order orientational order-disorder tran- 
sition [13] at T = T m ~ 260 K, the molecules actually unlock to attain seemingly 
free rotation at T > T m . At T m there is indeed an abrupt change in friction 
but the lowest friction coefficient is still of order /i ~ 0.15 0, 0], worse than with 
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traditional oil-based lubricants. Coffey and Krim [5] undertook a quartz crystal mi- 
crobalance study of one or two Ceo monolayers adsorbed on Ag(lll) or Cu(lll). 
There are no rotations in a C§q layer on Cu(lll), and only a slow change of molec- 
ular orientations for a Cgo layer on Ag(lll). In a bilayer, Cgo molecules in the top 
layer rotate freely at 300 K, as expected. Nonetheless a molecularly thin methanol 
film deposited over the Cgo bilayer failed to show either the expected low fric- 
tion, or any essential difference between these systems. Altogether, these results 
apparently doom the idea to replicate the bearing concept at the nanoscale. The 
questions why, and how this could be overcome are still quite open. 

First of all, even a small amount of charge transfer and/or bonding between Ceo 
and the metal substrate may hinder the rolling; therefore, chemically inactive, in- 
sulating substrates should be preferred to metals. A second handicap lies in the 
close packing of molecules. Balls in bearings are arranged so as to prevent con- 
tact, but rolling molecules in a monolayer are in close contact, hindering, even in 
the apparently freely rotating phase, their mutual rolling and jamming the same 
way mutually ingrained rolling cogwheels would. As was found in earlier simula- 
tions a way to avoid jamming is to lower dramatically the molecular coverage 
well below a monolayer - indicatively to coverages < 0.3 - to approach the limiting 
rolling friction of a single molecule. A third element, not investigated so far, is a 
possible role of molecule-substrate matching. Using a single molecule model as our 
starting point, we now find that this element is unexpectedly important. 

In our ordinary, macroscopic world, the main source of rolling friction is de- 
formation, since both substrate and roller elastically or plastically deform at the 
contact. The deformation energy is partly released and lost in the form of bulk 
frictional heat when the roller moves on PJ]. By desi gning the bulk so that dissi- 
pation is poor in the right frequency range, rolling friction can be made 10 2 to 10 3 
times lower than sliding friction, since the latter implies breaking and re-forming 
of slider-substrate bonds. As the roller size is decreased however, adhesion grows in 
importance, eventually becoming the main source of friction. To rotate a molecule, 
one has to break the molecule-substrate bonds one side of the molecule and cre- 
ate new bonds on the opposite side; hence, molecular rolling friction might not be 
much lower than sliding friction. 

Our goal is to understand what physics may yield the lowest friction coefficient 
attainable for molecular rolling and which system parameters might provide it. In 
the end, we find that nanobearings might indeed work as well as macroscopic ones, 
but one has to choose properly the macroscopic counterpart, which here turns out 
to be a perfect rack-and-pinion matching as in cogwheels. Because we are interested 
in general trends, we presently explore a minimal 2D model, which allows us to 
span a large number of parameters, and also provides an easier visualization of all 
processes. 

2. Model 

Consider two substrates with lubricant molecules in between, all elements made 
up of classical point particles (atoms). Atoms can generally move in the (x,y) 
plane, where x is the sliding direction and y is perpendicular to the substrates. 
The substrates, pressed together by a load force i 7 ) = N s fi, consist of rigid atomic 
chains of length N s and lattice constant R s , so that the system size in the slid- 
ing direction is L x = N S R S and the total mass of the substrate is N s m s (we 
use periodic boundary condition along x). The bottom rigid substrate is fixed at 
x = y = 0, the top one is free to move in both x and y, while driven along x 
through a spring of elastic constant k s moving with speed v s . The spring force 
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F, whose maximum value before the onset of sliding measures the static friction 
force F s , and whose average during smooth motion = (F) is the kinetic fric- 
tion force, is monitored during simulation (throughout the paper we will normalize 
forces per substrate atom / = F/Nc). Thus, our model is a 2D variant of a typical 
experimental setup in tribology [l5|. Between the substrates we place circular 
("spherical") lubricant molecules. As in Ref. 13 ] . each molecule consists of one 
central atom and of L atoms on circle of radius R m = Ru/2sm(n/L) so that their 
chord distance is Ru. All atoms interact via a 12-6 Lennard- Jones (LJ) potential 



^Lj(r) 



V„ 



(R aa '/r) 12 -2(R aa ,/rY 



where a, a' = s or I for the substrate 



or molecular atoms respectively. The latter are additionally coupled to the cen- 
tral atom by stiff springs of the elastic constant K m , I^tab^) = \K m {r — i? s tab) 2 , 
where the distance R sta ,h = R m + (12 V n /K m R m ) [(R u /R m f - (Ru/R m ) 12 ] is cho- 
sen so that the total potential VLj(r) + T4tab(0 is minimum at r = R m . With 
K m = 100 the resulting stiff molecular shape resisted destruction during the sim- 
ulations. Thus, the lubricant-lubricant interaction is described by the parameters 
Vu and Ru, while the lubricant-substrate interaction, by V s i and R s i (direct inter- 
action between the top and bottom substrates is omitted). We use dimensionless 
units, where m s = mi = 1, Ru = 1, and the energy parameters V aa > takes val- 
ues around V aa > ~ 1. Because a 2D model cannot reproduce even qualitatively 
the phonon spectrum of a 3D system, and because frictional kinetics is gener- 
ally diffusional rather than inertial, we use Langevin equations of motion with 
Gaussian random forces corresponding to temperature T, and a damping force 
fn,x = —mr](y)x — mrj(Y — y) (x — X), where x,y are the atomic coordinates 
and X, Y are the coordinates of the top substrate (the force f n>y is defined in the 
same way). The viscous damping coefficient is assumed to decrease with the dis- 
tance from the corresponding substrate, rj(y) = rjo [1 — tanh(y /yd)), where typically 
7]q = 1 and y d ~ 1. 



3. Results 
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Figure 1. Activation energy E a as a function of the substrate lattice constant R s for the rigid L = 6 
molecule, for /; = 0.5, V s i = 1/9, and R s i = Ru. Open symbols correspond to sliding, solid symbols to 

rolling. 
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Simulation results for single molecule friction were obtained from L = 5 (the 
simplest circular molecule) up to L = 13 and 14, which may be considered as a 2D 
version of fullerene. In fact in the 3D case, the surface area of the spherical molecule 
is s = 4irRm- With L3 = 60 atoms on the surface, this gives s L^Rf l7 or R m /Ru ss 
2.18. In 2D, the length of the circle is 2irR m ~ LRu, or L « 2irR m /Ru, which leads 
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a (min) b (max) c (min) d (max) e (min) 




Figure 2. (color online): Sliding adiabatic motion of the rigid L = 6 molecule for R s /Ru = 0.66 
(A</< < </<o, left panel) and rolling for R a /Ru = 1.29 (A<j> > <j>0> right panel). Other parameters as in 
Fig.[TJ Lower panels: X-dependence of the potential energy V(X), and the angle <f>(X)/<t>0- Top panel: 
configurations as the molecule moves from one minimum of V(X) to the next. 




to L ~ 13.7 for the same ratio R m /Ru as for 3D fullerenes. We first consider a rigid 
molecule, Vu = oo and K m = oo, fix X of the top substrate and seek the potential 
energy minimum V by varying its y-coordinate Y, the molecular center (x c , y c ) and 
the rotation angle 4>. The X dependence of V, (x c ,y c ), and (f> defines the adiabatic 
trajectory, which describes the joint substrate and lubricant motion when infinitely 
slow. We extract the activation energy E a = max [V(X)] — min[V(X)], and the 
static friction force magnitude, approximated as f s = max [dV(X)/dX] (f s ~ E a 
in our units). 

Figures [T]and [2] show results for the L = 6 molecule when, to simplify further, R s i 
is fixed to R s i = Ru. The energy V(X) is periodic with R s (or a multiplier of R s ). 
The molecular angle eft varies by A(j) as the potential energy V(X) changes from 
minimum to maximum. Because <ft(X) must be continuous, the motion corresponds 
to sliding if Acp < 4>o = 2tt/L, while if Acp > 0o the molecule must rotate when 
it moves. As Fig. Q] shows, for R s < Ru the motion corresponds to sliding, i.e., 
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r ♦ data from Fig. 3 

fit by Eq(l) with ^=0.31 (evenl) 

1 4 I" fit by Eq(l) with ^=0.17 (odd L) 




Figure 4. (color online): Position R s of the main minimum of E a (R s ), extracted from Fig. [3] for 
increasing molecular size L. Curves are fits to the cogwheel model Eq. JTJ. The asymptotic limit of 1 is 
still relatively far for reasonable molecular radii. 




Figure 5. (color online): Static friction force /„ as a function of the substrate lattice constant R B /Rn for 
(a) L = 14 and (b) L = 13 for different system parameters: (i) /j = 0.5 and V 3 i = 1/9 (solid curve and 
circles and stars), [it) fi = 0.5 and V 3 i = 1/3 (down triangles and red dotted curve and stars), and 
(in) fi = 0.1 and V a i = 1/9 (up triangles and blue dotted curve). Other parameters are: R a i = R B , 

K m = 100 and V H = 1. 



the molecule is shifted as a whole, slightly oscillating during motion (see Fig. [21 
left panel). The sliding activation energy has maxima at Ru = nR s (where n is an 
integer) and minima at Ru = (n — 1/2)R S . On the other hand, for Ru < R s the 
motion corresponds to rolling (Fig. [21 right panel). Here E a (R s ) has sharp minima 
at unanticipated values of R s /Ru (e.g., for R s /Ru « 1.29 in Fig. [T]). 

Varying R s in Fig. [Hat fixed value of R s [ for the lubricant-substrate interaction, 
we find that for R s < Ru rolling replaces sliding already around R s /Ru = 2/3. 
Preference for rolling over sliding increases for increasing load /; and for decreasing 
interaction strength V s i (choosing alternatively R s i = R s rolling is even more 
prevalent than for fixed R s i)- When sliding wins over rolling, it generally provides 
a lower activation energy. Recalling that 4>o = 2ir/L, the region of parameters for 
rolling should increase with L — a rounder wheel rolls better! The R s dependence 
of E a (R s ) for increasing size L (Fig. [3]) shows rolling for all R s and for all L > 5, 
except for L = 6 which shows both rolling and sliding (see open symbols in Fig.[3h). 
As R s varies, the value of E a oscillates by more than two orders of magnitude for 
even L and more than three for odd L, with deep sharp minima separated by 
broad maxima. Clearly, by suitably choosing R s /Ru a very strong decrease of 
rolling friction is attainable. 

The unexpected minima of E a {R s ) can be explained by simple engineering — a 
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"rack-and-pinion" model. Consider the molecule as a cogwheel (the pinion) with L 
cogs, primitive radius R m and external radius R* = R m + h, where h oc R s i. The 
chord distance between nearest cogs is = 2R* sm(ir/L). Best rolling conditions 

are expected when R^ matches the substrate potential period R s , i.e., for Rs = R* u 

(2) (3) 

and its fractions, Rl 7> = R* u /2, R\ = R^/3, etc. The main minimum of E a {R s ) is 
expected at 

R^/Ru = 1 + (2h/R u ) sm(n/L). (1) 

As shown in Fig. HI the cogwheel model ([1]) with h = (3R s i, where (5 is a parameter, 
fits extremely well the shift of minimum position with molecular size L. Moreover, 
it can explain its variation with load (the radius R* and therefore h decreases as 
the load grows) as well as with the lubricant-substrate interaction V s i (R* and h 
decreases with V s i) [lfj]. It also accounts for the even-odd effect since odd L involves 
perfect ingraining for one substrate at a time, justifying why roughly double values 
of /3 are needed for even relative to odd L. 

These insights for rigid molecules are confirmed by the static friction force ob- 
tained from simulation with deformable molecules (Fig. [5]). The friction coefficient 
Ms = fs/fl ranges from fi s ~ 0.1 at R s /Ru 0.7 to fj, s ~ 0.01 or even fj, s ~ 0.001 
at R s /Ru ~ 1.1. The results are also robust to changes of model parameters. For 
example, Fig. [5] compares the dependences f s (R s ) for two values of the amplitude 
of lubricant-substrate interaction, V s i = 1/9 and 1/3, and for two values of the 
load, fi = 0.5 and 0.1. We found both static f s and kinetic friction to increase 
approximately linearly with load, f S; k ~ fos,ok +fJ> s ,kfl- Visualization of MD trajec- 
tories shows that for R s /Ru = 0.7, where friction is high, rolling rotation is in fact 
discontinuous, and accompanied by a molecular shift /sliding — much as cogwheels 
with excessive clearance would do — which dissipates mechanical energy into vi- 
brations, whereas for R s /Ru = 1.1, where friction is low, the motion is a smooth 
rotation, corresponding to optimal rack-and-pinion coupling. 

Simulations further showed that this scenario remains valid at nonzero tempera- 
ture T. As T increases, both static and kinetic friction forces are found to decrease, 
the stick-slip changing to creep and finally to smooth motion at high tempera- 
tures. We also found transitions from stick-slip to smooth rolling for increasing 
velocity. The cogwheel effect remains, and for example calculated static friction for 
R s /Ru = 0.7 and R s /Ru = 1.1 still differ by a factor of ten or more. The critical 
velocity v c of the transition from stick-slip to smooth rolling also differs by a factor 



of about four in the two cases. In all cases we find fk <C f s ly |. 

The present approach to the single rolling molecule can be extended to a finite 
coverage of lubricant molecules. For M molecules the individual molecular load is 
fli = fi/M so that the friction per molecule f s i,ki a ^ so decreases. However, the 
total friction force f s ^ = Mf s \}.\ increases, so that the combined effect is a slow 
increase of the total friction with M. In a real situation, coalescence may lead to 
jamming, with molecules blocking their mutual rotation. In our model, jamming 
starts already at 9m ~ 0.1 and completely destroys rolling at 6m > 0.3 13|,lla] (here 
6m = M/M\ is the coverage, with M\ the number of molecules in a monolayer). 



4. Conclusion 

Our 2D modeling leads to overall conclusions of considerable novelty. Rolling spher- 
ical lubricant molecules can indeed provide better tribological parameters than 
sliding atomic lubricants. The effect may be as large as in macroscopic friction, 
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where rolling reduces friction by a factor of 10 2 — 10 3 , however only for low concen- 
tration of lubricant molecules, and for specially chosen values of the ratio R s /Ru, 
e.g., for R s /Ru ~ 1.1, corresponding to perfect rack-and-pinion matching. 

While of course the matching condition will be more delicate and difficult to 
realize in 3D than in 2D, the concept can surely be pursued experimentally, for 
different spherical molecules, different substrates, and different coverages. Inert 
nonmetal surface (such as perhaps self-assembled monolayers) could represent a 
better choice of substrate than metals for fullerene deposition; and lower coverages 
should be preferred to complete monolayers. Because 3D rolling has an azimuthal 
degree of freedom, the novel cogwheel effect found here will be direction dependent, 
and rolling friction will exhibit anisotropy depending on direction. At fixed rolling 
direction, increasing load could offer the simplest tool to change the rack-and- 
pinion matching through a pressure- induced decrease of h = (3R s i- In this case, 
the effect would show up as a strikingly nonmonotonic (non-Amontons) behavior 
of friction with load. 
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